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, The partition function of a two dimensional Abelian gauge model reproducing magnetic vortices 

^\ ■ is discussed in the harmonic approximation. Classical solutions exhibit conformal invariance, that 

0^ ' is broken by statistical fluctuations, apart from an exceptional case. The corresponding "anomaly" 

has been evaluated. Zero modes of the thermal fluctuation operator have been carefully discussed. 
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I. INTRODUCTION 

In the recent years, a great interest has been devoted to the study of magnetic vortices, which are intimately 
connected to the study of superconductivity . A magnetic vortex can be represented as an infinitely long magnetic 
. ' flux tube in three spatial dimensions or, equivalently, as a two dimensional localised magnetic field source. The first 
vortex solutions were discussed in 1957 by Abrikosov Q in the framework of the Ginzburg-Landau model, whereas 
5_i ' Nielsen and Olesen [|| in 1973 were the first to recognize that the same type of vortices are also present in a classical 
relativistic model, i.e. the Abelian Higgs model. Several years later Hong, Kim and Pac 0] and, independently, 
Jackiw and Weinberg ||, discussed vortex solutions in the same (planar) model, but with dynamics for the gauge 
fields governed by a Chern-Simons (CS) term. In a CS system, the magnetic field is proportional to the charge 
density, so that any excitation carrying magnetic flux is necessarily charged, contrary to the vortices in the Abelian 
Higgs model, that are electrically neutral. 

Unfortunately, in both cases the classical solutions are known either asymptotically or by series with complicated re- 
cursion formulas ^ . Consequently, quantization or thermal fluctuations of the system around the classical background 
are very difficult to study. 

More recently, a different type of magnetic vortex was proposed Q as a solution of 2 dimensional euclidean scalar 
electrodynamics with topological coupling. In such (Liouville) vortices, the magnetic field satisfies the Liouville 
equation, all of whose regular solutions can be easily expressed in terms of an arbitrary analytic function. In spite 
of the seeming resemblance of this model with the Abelian Higgs system, Liouville vortices are very different from 
the ones in the Ginzburg Landau theory. In the Ginzburg Landau case, the potential of the Higgs field leads to a 
spontaneously symmetry breaking. Due to this fact the classical vortex solutions, besides the expected long range 
tail, have an exponentially decreasing component, and the coefficient of the quadratic term in the Higgs potential is 
related to the characteristic length of the exponential behavior. On the contrary, in Liouville vortices the potential 
of the scalar field is that of a pure \4>\ theory, without quadratic-mass term. In addition, there is also a non-minimal 
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(topological) interaction that couples the matter density directly to the magnetic field. As a consequence of these two 
facts, Liouville vortices exhibit conformal invariance, and their asymptotic behavior is always inverse power-like. 

The model leading to Liouville vortices shares important properties with two other models: the Jackiw and Pi model 
H and the non linear sigma model (NLSM) with 0(3) local symmetry M. Concerning the first, the profile of the 
magnetic field of these Liouville vortices is identical to that discussed by Jackiw and Pi as the static soliton solution 
of the gauged non linear Schroedinger equation on the plane, in strong analogy with the Abelian Higgs model, that 
is a classical field theory whose equations of motion coincides with the non linear Schroedinger equation governing 
the Ginzburg-Landau theory for type II superconductors. Concerning the NLSM, in ref. it was shown that all the 
solutions of the above mentioned 2 dimensional euclidean scalar electrodynamics can be obtained from the solutions 
of NLSM with local symmetry. In addition, there are other important features shared with the euclidean NLSM in 2 
dimensions Q that will be discussed below. 

Statistical mechanics is a natural framework for topologically non trivial euclidean theories and, on the other hand, 
many important physical features that the model in ref. should hopefully exhibit, are strictly related to statistical 
mechanics. As an example, if these Liouville vortices are really related to superconductivity, a phase transition of the 
system should occur at some critical temperature. 

In euclidean models, with positive defined actions S and admitting topologically non trivial solutions, the partition 
function can be defined as the path integral over the configuration space of the Boltzmann factor e~ s In this 
context, the action plays the role of potential energy and the free energy is usually interpreted as interaction energy 
due to "thermal" fluctuations. However, in these cases, the definition of temperature in not always straightforward. 
For instance, in the NLSM it is common to consider the classical action already implemented by an overall factor 
1/e 2 , i.e. 

S = 2? / d2xd » Nad n N<1 > NaNa = 1 > (!) 

and such a factor is eventually interpreted as "thermal bath" (3 = 1/e 2 . 

As we shall see, in our model it will be possible to rescale the fields in such a way that the action depends on 
the U(l) coupling e only through an overall factor 1/e 2 , just like in the NLSM above, allowing a temperature-like 
interpretation of the coupling constant. 

The great advantage of such Liouville vortices is that the classical solutions are given by relatively simple expressions, 
so that the problem of the thermal fluctuations of the classical solutions can be faced. To provide a first insight in 
this direction is the aim of the present paper. 

In Sect. II we shall briefly review the model. The solutions of this model are conformally invariant and can be 
classified according to their vorticity, which is proportional to the topological invariant of the model and is nothing 
but the degree of the arbitrary analytic function tu(z) upon which the solutions depend. 

In Sect. Ill we begin the study of the thermal fluctuations of the model. After noticing some remarkable similarities 
with the euclidean NLSM in 2D, we consider the partition function in the so called harmonic approximation: fields 
are parametrized as classical solutions plus small thermal fluctuations (small coupling regime) , and then the action in 
the Boltzmann factor is replaced by its expansion around the classical solutions up to the second order in the thermal 
fluctuations. All the zero modes of the corresponding fluctuation operator can be explicitly evaluated, for any choice 
of the arbitrary analytic function lo(z) characterising the background. 

If the vorticity J\f of the background is different from the minimal one (i.e. Af 1), the conformal symmetry of 
the classical solutions is broken by thermal fluctuations, and in Sect. IV the corresponding "conformal anomaly" is 
explicitly calculated. 

In Section V we discuss the special case M = 1 and we comment on future developments. 



II. THE CLASSICAL MODEL 



Let us consider the 2-dimensional euclidean action of a |</>| 4 complex scalar field minimally coupled to a Maxwell 
jauge potential with an additional topological coupling 



d 2 x 
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-iiijF^v + ^ ( 



(2) 



where D^cfi = d^ifi — ieA^cfi is the U(\) covariant derivative and e is the Abelian coupling constant that in 2D has 
dimensions 1 in mass units. The topological coupling ecf>* (fie ^ v F^ v is consistent with both C7(l) gauge invariance and 
ISO(2) invariance of the action (|2|). Actually, any general coupling of the type ef((j>*<fi)eijFij and any general scalar 
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potential V((/>*(/>), f and V being arbitrary functionals of the scalar density 4>*4>, would be also allowed, due to the 
fact that in 2D the scalar field <f> is dimensionless. Such a general case has been considered in ref. jjj and, with a 
suitable choice of the two functionals, it is possible to obtain classical solutions of the generalised version of the action 
(|J) satisfying a wide variety of 2 dimensional conformal equations. Here we shall restrict to the most interesting case 
given by (§). 

Up to a total derivative, S can be rewritten as 



S= / cfx£ = 



{D x ±iD y ), 



> , (3) 



B = F\2 being the magnetic field. As a consequence, (|^) is extremized by field configurations satisfying the following 
self-duality (anti self-duality) conditions 

(Dp ±ie llv D v ) = 0, 

B = ±ecj)*(t). (4) 

It is easy to verify that any configuration satisfying (jl]) also solves the classical Euler-Lagrange equations. Combining 
together eqs. (||), the scalar density 0*0 satisfies the Liouville equation 

A In (0*0) = — 2e 2 ((j>*(j>) , (5) 

all of whose regular, positive definite solutions are given by 

4 M 2 



[l + \w(z)\ 



(6) 



In eq. (Q), w(z) is an arbitrary meromorphic function and u)'(z) — dui/dz. Explicit solutions of the self-duality 
conditions (|4|) in the covariant gauge d^A^ = are 



e A^ = -e^ v d v In [l + | oj (z) | ; 
,_ 2a/(z) 
^ 1 + M*)l a ' 



(7) 



The above field configuration describes a magnetic vortex. As a matter of fact, for any choice of the arbitrary 
meromorphic function lu(z) the magnetic field is always localised. The magnetic flux &(B) is quantized, as the 
integral of the magnetic field over the whole space is proportional to the degree of the analytic function oj, i.e. the 
number of solutions Zi = Zi(u>) of the equation lu = lu(z), each multiplied by the appropriate multcplicity hi, namely 



The integer Af is usually denoted as vorticity, and is the topological invariant associated to the classical configuration. 
For the self-dual configurations we are considering, Af > 0. A parity transformation on the solutions maps self-dual 
into anti self-dual configurations and the vorticity changes sign. 

In 2 dimensions the natural flux units are $o = 27r/Q, Q being the electric charge. Consequently, if we identify e 
with the electric charge Q, the flux is an even multiple of <I>o. Alternatively, eq. (g) gives a magnetic flux which is 
an arbitrary integer Af in terms of the natural flux units $o > provided the Abelian coupling e is twice the electric 
charge Q , suggesting the idea that the matter field we are considering should be somehow related to an electron-pair 
condensate (Cooper's pair). 

In the remaining part of this Section, we shall discuss the classical symmetries of the solutions (^), as they will play 
a crucial role in the thermal fluctuations the model (^|). Besides the obvious gauge and 'Poincare' ISO(2) symmetries, 
solutions (|7|) possess also conformal invariance. A conformal transformation 

z = x + iy -> p(z) = x(x,y) + iy(x,y) , (9) 

where p is an arbitrary analytic function, connects different solutions of the action (Q). Under the action of the 
conformal redefinition (^|) , gauge and matter fields transform as 
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cj)(z, z) -> 4>(z, z) = j^4>{p, p) 
^(r)^A>) = g^(r) 



(10) 



so that the matter density </>*</>, the magnetic field and the self-dual derivatives transform as densities with weight 
J = dct(9^ / r dx v ) = \p'{z)\ 2 , namely 



(j>*(/) ~ > (j)*4>(z, z) = J(j)*(j) 



[(I> M + ie^DJ) 0] (r) 



Bv + iea V D v 



(p,p) 

(r) = J [(£>„ + ie^M 0] (f ) 



-> B(z,z) = JB 



(p,p) 



(11) 



^From eqs. ( |ll| ) it follows that both conditions (|4|) are preserved under a conformal transformation. Conformal 
symmetry of the solutions can be easily understood in terms of the energy momentum tensor T M „ . The independent 
components of in 2 dimensions can be chosen to be its trace T M/J and = T\\ — T22 i 2zTi2 . In our specific 
model p|) such components are given by 



T W - o ( B 



(B 



T ± = 2 (£>!</.* ± i£> 2 0*) (-Di0 ± iB 2 <; 



(12) 



and we see that the self-duality conditions ([|), that solve the equations of motion, also make vanishing the whole 
energy momentum tensor, and in particular its trace. 

Nonetheless, the action S is not conformally invariant. At first sight, conformal symmetry of the solutions could 
seem at odd with the fact that the action explicitly depends on a dimensional parameter (the Abelian coupling e), and 
with the fact that the action is not conformally invariant. However, to preserve the stationary points of the action, 
the Lagrangian density transforms covariantly under conformal transformations, 



d 2 x £(x) 



(fxj£(x) 



(13) 



The reason of such a nice transformation property relies in the fact that it is always possible to rescale the fields in 
such a way that the action (||) does not contain any mass scale, apart from an overall multiplying factor, just like in 
the NLSM (0). As a matter of fact, if we rescale the fields, collectively denoted by as $ = e<i>, the classical action 
(B) written in terms of $ becomes 



(14) 



and, classically, all the dependence on dimensional parameters in the action (0) can be ruled out. 



III. STATISTICAL FLUCTUATIONS 



The partition function of an euclidean system admitting topologically non trivial solutions is defined as the path 
integral over the configuration space of the Boltzmann factor e~ s |9|]. 

There are some important analogies between the model we have described in the previous Section and the NLSM 
in eq. ([!]) that are definitely worth mentioning, also in view of the fact that for the NLSM the partition function can 
be estimated in the small coupling regime ||. 

In both the models, all the classical solutions are expressed in terms of an arbitrary analytic function w(z), and 
the topological invariant expressed in terms of io{z) is identical for the two models. In addition, due to eq. ([il|), 
also in our model all the dependence on the coupling constant can be factorized in an overall factor 1/e 2 , allowing 
a temperature-like interpretation of the coupling constant in the evaluation of the partition function, just like in the 
NLSM. 

Finally, the most remarkable analogy: following in order to exploit the renormalization group of the system (Q), 
we decompose the N a variables of the NLSM according to 

7V a (a ; ) = (l-|^| 2 ) 1 / 2 iV a (^) + ^ , (15) 
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where a = 1,2,3 and i — 1,2. Nfi is the so called "slowly varying" vector and ef are orthogonal to it, whereas 
represents the "fast fluctuations". As a consequence of iVge" = and e"e° = <£y, we have 



= B^ef 



d e a 



B^NS . 



(16) 



The variables and have to be considered as auxiliary variables characterizing the slowly varying fields. Substi- 
tuting eqs. (15, [Lq) in cq. ffl) and selecting terms up to the second order in <p l one gets 



d 2 x 



(17) 



where <fi — (0 1 +*0 2 )/ V% and F Mt/ is the "field strength" of the auxiliary variable A^. Notice that eq. ([Tt]) is exactly of 
the same type of the quadratic term in <j) of the action (g); in particular, the same topological coupling is reproduced. 
However, its coefficient is one half of the one in eq. (||), and therefore eq. ( |l7| ) is nothing but the first term in the 
r.h.s. of eq. (||), up to an inessential total derivative. Consequently, extremizing the action S 1 - 11 ^ is equivalent to 

(and the second set of 
has to be considered as 



impose the first set of self-duality conditions (|4|). The remaining term of the action in 
self-duality conditions (Q)) provides the "dynamics" for the field A^ that, instead, in eq. (|17 
a background field. 

Thus, a first very rough estimate of the partition function of the model (||) could be obtained in the following way. 
In the path integral, one could integrate only over matter fluctuations (p — — </> c ' induced by the first term in eq. (|^), 
keeping the magnetic field fixed and equal to its classical value B = (f>* cl (f) cl . Then, one immediately gets Polyakov's 
results and the partition function turns out to be that of a Coulomb gas in its plasma phase, with Debye screening, 
and therefore a mass gap. 

On the one hand, thefact that the original system is equivalent to a set of massive fermions is definitely a positive 
result, that confirms previous conjectures 0. On the other hand, a physical interpretation of this picture in terms 
of vortices is still obscure and, more important, it is not clear to what extent the approximation of considering only 
matter fluctuations is reliable. 

To proceed, one should consider fluctuations of the whole set of self-interacting fields (matter and gauge fields). 
Clearly, such a problem is much more complicated and an explicit evaluation of the partition function becomes a 
formidable problem. However, even without an explicit knowledge of the partition function, something can be done 
and interesting results can be obtained. 



A. Harmonic approximation 



In the evaluation of the partition function of our model, the most appropriate approximation is the so called 
harmonic approximation, where the action in the Boltzmann factor is expanded up to the second order around the 
classical solutions. This approximation is frequently used also in quantum mechanics (stationary phase approximation) 
where it leads to the one loop effective potential. 

Clearly, in the harminic approximation, it is assumed that trajectories deviating significantly from the classical 
solutions have a negligible weight in the path integral. In our case, due to eq. (|l4|), this assumption is certainly true 
in the small coupling regime. Thus, we shall consider the case f3 — 1/e 2 ^> 1. 

Let us decompose the fields $ = {Ai, A2, (f>i, $2}, with </> = (l/V2)(</>i + ifa), as the sum of the classical solutions 
$ d = {Af , AS? , <t>f , 
expansion of the action 

two terms of this expansion vanish, due to the Euler-Lagrange equations and to the fact that (|2|) vanishes on the 
classical solutions. We have therefore 



plus "thermal" fluctuations r\ = {ai , a% , <pi , ^2} and let us consider the formal 
I) around the classical solutions up to the second order in the field fluctuations. The first 









4/ 



5 2 S 



S 

A ^(x)^-(y) 
Fluctuations 77 are required to be sufficiently regular and normalizable, 

||ry|| 2 = / d 2 x [a M a M + tpwi] < 00 



(18) 



(19) 
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Applying eq. ( |l8| ) to our model, after straightforward calculations, we get 



S|$J ~- / d 2 xd 2 yr h (x)5 2 (x-y)M ij (y)r ]j (y) 




= - / d 2 y { a M [(-A + e 2 (0 2 + 2 )) 5^ + fl^ft,] a v + 

+a M [-2 e 0i e^dy - e e^d^x - 2 e 5^02 + 2 e 2 A^4>i] ipi + 

+cpi [2 e 0i e^c^ - e e^d^i - 2 e <9 P 2 + 2 e 2 A^ 0i - 2 e fad^] a p + 

4a M [-2 e 02 £/uA - e e^d^i + 2 e 9^01 + 2 e 2 A M 2 ] ^2 + 

-H^2 [2 e 02 e M iy9r/ - e e^dvfc + 2 e <9 p 0i + 2 e 2 A M 2 + 2 e0i<9 M ] a M + 

-Hpi -A + i e 2 (0 2 -0^) + e 2 A 2 (^i + (^i [e 2 0i02-2eA Ai a A1 ] <^ 2 + 
-H^2 [e 2 0i0 2 +2eA fl d fl ] (pi + tp 2 - A - - e 2 (0 2 - 0^) + e 2 ^ 2 <^ 2 }, 



(20) 



where A 2 = A^A^ and we omitted, for brevity, the superscript cZ on the classical fields in the square brackets. In 
order to factorize the 5 2 (x — y) term in the first equality of ( po| ) , some integrations by parts have been performed 
and, in so doing, the fluctuation operator Mij does not look manifestly self-adjoint. This is not a problem, as the 
normalizability condition (|l^) always permits to write as in ( |20| ) . 

Evidently, in the harmonic approximation the evaluation of the partition function is equivalent to the calculation 
of the determinant of Mij and, in turn, the solutions of the eigenvalue problem 



would completely solve the problem. 

Clearly, this approach is practically impossible, due to the complicated form of the operator My. Nonetheless, even 
without the explicit knowledge of the determinant, quite often it is possible to calculate expectation values. This is 
the case of the conformal anomaly, that we shall evaluate in Sect. V with the help of the zeta-function regularization 
technique. 



^From eq. S is manifestly positive definite and, on the classical solutions, achieves its vanishing minima. 
Consequently, has no negative eigenvalues, and the lowest eigenvalues are the zero modes. In our formulation, 
there are many operators My, depending on the choice of the classical background that, in turn, is completely specified 
by the choice of the arbitrary meromorphic function lo(z). For a given choice of ui one can associate a topological 
invariant J\f (vorticity) to the classical solutions. The topological invariant classifies classical vortex solutions into 
distinct inequivalent classes. Obviously, in each class there exist infinite analytic functions leading to the same 
vorticity. 

Zero modes will strongly depend on the specific choice of the arbitrary analytic function u> characterizing the 
classical background. Nevertheless, the number of normalizable zero modes of the operator My will depend only on 
the homotopy class to which the function uj belongs. 

Zero modes satisfy a first order equation which is simpler than eq. (pTI) with A< = 0, and therefore is worth 
mentioning. Such an equation is a consequence of the peculiar form (0) of the action. Due to the self-duality 
conditions ([|), the action vanishes when evaluated on the classical background (R). On the other hand, from eq. (|^) 
the action is written as the sum of squares, so that the requirement that the action vanishes up to the second order 
in the "thermal" expansion $ = $ c ' + n (zero- modes), is equivalent to the requirement that zero modes n solve the 
self-duality conditions (0) expanded up to the first order, i.e. 



(21) 



B. Zero modes of the thermal fluctuation operator 



(D C J; + ie^ v D c j) p = ie (a M + ie^a^) 0' 



cl 




(22) 



G 



where D C J — — ieA c J , b = e^ v d^a v . It can be verified by direct inspection that eqs. (^2|) are indeed equivalent 
to eq. ( [HI) with Xi = 0. Moreover, eqs. (22) can be decoupled and rewritten as a Schroedinger type problem. Taking 
eqs. (||) and (||) into account, it is not difficult to check that the magnetic field b of the zero modes has to satisfy the 
equation 



b 
B 3 



= -2eb (23) 



that, written in terms of bj B cl — becomes the zero-energy Schroedinger equation for a unit mass wavefunction ip 
moving in the classical potential V = —eB cl , i.e. 

- ^Aip + Vip = , V = -eB cl . (24) 

Consequently, by finding the zero modes of our system one gets, as a by-product, the zero energy solutions of the 
above Schroedinger equation. 

Zero modes of the fluctuation operator M^- are associated to the continuous symmetries of the classical solutions. 
Such symmetries are: translations, rotations, gauge transformations, conformal transformations and variations of the 
arbitrary parameters upon which the classical solutions may depend. Practically, once a gauge choice has been picked, 
all such symmetries are included in conformal transformations. 

It can be shown that, starting from a classical solution $ c/ , the zero mode of the operator My associated to 
the continuous symmetry E with infinitesimal parameter a is given by the variation of $ ci under the action of the 
continuous symmetry E , i.e. if $ c/ — * <& d + fe<I> d , then 77 s = S^ cL /Sa is the zero mode of associated to the 
symmetry E. The proof can be easily obtained by performing a transformation S on the Euler - Lagrange equations. 

The same procedure can be also generalised to local symmetries, like the gauge and the conformal ones. In this 
case, zero modes are obtained by performing an infinitesimal transformation on the classical fields: the obtained result 
is a zero mode also when the function specifying the continuous transformation is no longer infinitesimal. Several 
examples will be provided below. 



C. Explicit evaluation of zero modes 

We begin by evaluating the zero mode associated to the gauge symmetry. By performing an (infinitesimal) gauge 
transformation on the classical solutions we get 

r, G = (-^a; -a<$; a<j>^ (25) 

It can be verified that ( |25| ) is a solution of the coupled equations (^2|), and therefore a zero mode of My for any 
arbitrary function a (x) , not necessarily infinitesimal. However, the simultaneous requirement of normalizability (M) 
and the gauge condition d^A^ = forces a to be a constant. Introducing for later convenience complex notation 
a = a\ + i ei2 for the fluctuations of the gauge potential, the only normalizable zero mode associated to residual gauge 
symmetry is then 

77 G = (^, a ) = (ia0 ci ,0) . (26) 

We now consider the remaining symmetry transformations of the classical solutions. Such symmetries are all 
included in the conformal one, since traslations, rotations and variations of the parameters upon which the classical 
solutions may depend can be always seen as particular cases of conformal transformations. Actually, even the zero 
mode ( p6| ) arises from a particular conformal transformation. 

Since the classical solutions are completely determined in terms of the arbitrary analytic function lo(z) (see eq. 
we can write a general expression for the zero mode associated to a given symmetry E in terms of Auj and u>, 
Aoj being the variation of uj under the action of E. By direct calculation we obtain the following general form of the 
zero modes 77 = (ip, a): 



e a = 2id 



2(Aw)' 2w'(wAw + wAw) 
1 + M 2 (1 + M 2 ) 2 



6( f= 1 I I. .19 T, I L. m2 ' ( 2? ) 



l + M 
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Now it is immediate to check that the zero mode ( |2q ) is a particular case of (|27]) , with Aw — iuj . In turn, such a Aw 
can be always obtained through a conformal transformation, so that hereafter the zero mode ( ^7|) will be classified 
among the conformal ones. 

The modes (^) are already in the Lorentz gauge (da + da — in complex notation). Among them, we have to 
select only the normalizable ones by imposing eq. (|l9|). In turn, normalizability condition ( |l9| ) can also be written 
in terms of w and Aw. After straightforward calculations, one can see that such a requirement is equivalent to the 
convergence condition of the two following integrals 



h = 
h = 



d 2 x 



d 2 x 



(Aw)' 



< oo 



(28) 



w(Aw)' - w'Aw 



l + M 



< oo 



Equations ( |27| ) and fl28| ) define the zero modes of the fluctuation operator associated to the continuous symmetry 
E. It should be noticed that (^7j) are indeed eigenvectors of the operator with vanishing eigenvalue even without 
specifying neither the form of w, nor its variation Aw. On the contrary, the normalizability criterion (|2^) will depend 
on the particular choice of arbitrary function w as well as on the variation Aw associated to the zero mode. As a 
consequence, to continue, we have to provide some specific examples. 

Let us fix the topological number of the classical solution to be AT . Clearly, there are infinite functions w with such 
a degree. Here, we shall consider two limiting examples, in such a way that the evaluation of zero modes in all the 
other possible choices of w will be easily understood as intermediate between these two cases. 

The first example is the totally degenerate case, where all the vorticity Af is carried by a single vortex that, for 
convenience, will be located at the origin. Then, 



w(z) = - 
Vzo 



(29) 



where zq is a scale introduced to render w dimensionless, as required. The choice ( |29| ) corresponds to the radially 
symmetric classical solutions 



iA cl = 2iM 



'o 



eB c 



AN 2 



l+(^ 

r 



2 A" 



2 AT 



(30) 



Af 



where r = \z\ and ro = \zo\- In this totally degenerate case, all the zero modes are easily obtained by performing a 
conformal transformations of the type z —* z + x( z ) on the classical solutions. The corresponding variation on w is 
thus 



Aw = uj'(z) X (z) = — (-) x(z). 



(31) 



We have now to investigate on the form of the function x( z )- Since the classical solutions (pG) vanish only at the origin 
and at infinity, one can easily realize that the only possible singularities of the function x{z) can be at the origin and 
at infinity: a singularity of x i n an y other point C ^ would necessarily render the corresponding zero mode singular 
and not normalizable in z = £ . On the contrary, singularities of x a t the origin are admissible, provided that when 
inserted in (|27]) and (^8|) they cancel against the corresponding zeros of the classical solutions and of the integration 
measure. Thus, in this totally degenerate case, it is not restrictive to consider functions x as pure powers of the type 
Xn — Qn(z j ro) n , with g n complex coefficients, and count the number of independent functions Xn that render the zero 
modes normalizable. For a given x.n{z) and with w(z) of the form (E9), the zero modes (ETT) take the form 



2iM {r/r f M 
r 2 1 + (r/r ) 2 ^ 



ro 
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l + (r/r ) 2 ^" V^W ^ 1 + (r/r ) 



n — 1 » j- / \ n — 1 

Z 



Af \ (z\ Af , . 



27V 1 



^ r 2 l + (r/r )^ U ' 



.z 




/V(r/r ) 2AA _ fz^ n - X 



(32) 



Considering the asymptotic behaviour of the zero modes (|32j) at the origin and at infinity, one can easily see that eq. 
( |32| ) define normalizable zero modes as long asl — Af < n < 1. Thus, there are Af + 1 admissible values of n. Since 
the coefficients g n are complex, two linearly independent zero modes |l2] ] correspond to each value of n , and eq. ( [32] ) 
with 1 — Af < n < 1 defines the 27V + 2 zero modes of the totally degenerate case. 

Notice that zero modes associated to IS 0(2) symmetry and to a variation of the parameter z in (|2^) are included 
in eq. ( j32|) with n = (traslations) and n = 1 (rotations if g\ is purely imaginary, and variation of the parameter zq 
if g\ is purely real). In addition, from eq. (|3l|), the gauge zero mode ( p6| ) is also a particular case of eq. (^) with 
n = 1. 

As a second limiting example we shall consider the totally non degenerate case, i.e. the choice of oj that, for a given 
value Af of vorticity, depends on as many free parameters as possible. Such an ui(z) is for instance, 

i— 1 v 7 

with a, 7^ aj, 6, ^ bj if i =/= j and a, 7^ bj for any i, j. It depends on 4Af + 2 real parameters (the constants ojq, 
{a,i} and {bi} are complex) and it corresponds to a classical TV-vortex configuration, where each vortex carries the 
minimum vorticity. The explicit location of the vortices is irrelevant to our purposes, but in general it will be a 
function of all the parameters juj . 
Let us define 

A=Y[(z- ai ), B = JJ (* - h ) , ( 34 ) 

i=l i=l 

in such a way that w(z) and its variation under a conformal transformation z — > 2; + x( z ) can be rewritten as 

w O0=^o^, Aw = w'(z)x(2)=wo^X! f^-^--^-^-Jx(^) ■ (35) 

i=l ^ 1 1 ' 

We have to investigate on the possible form of the arbitrary function xi z ) associated to the conformal variation 
Aw that makes the integrals ( |28| ) convergent. We first notice that, if x( z ) is regular in z — dj, then also (Aw)' is 
regular in a;. Thus, integrability of I± and I2 in z = only requires x(z) to be regular at 2 = a,. In z ~ b{, (Aw)' 
behaves like x(bi)/(z — bi) 3 . On the other hand, the quantity (1 + |w| 2 )~ 2 goes to zero for z ~ 6; like \z — bi\ 4 . Thus, 
in order to render convergent I\, x( z ) must have at least simple zeros in z = bi, and therefore it has to be of the type 
x( z ) = B f{z), with / regular in bi and a%- Finally, integrability of I\ at infinity requires f(z) ~ l/z q with q > Af — 1. 
There is only one possibility to have such a behaviour without introducing extra singularities in Aw, that is when the 
poles of f(z) exactly cancel with the zeros of ui'(z). Consequently, f(z) and x{ z ) have to be of the form 

fM z ) = 7~ — » XM{z) = BfAr(z) , (36) 



where Pj^-(z) is an arbitrary polynomial of degree A/". It is not difficult to verify that for such a choice of x( z ) & l so 
I2 is convergent. Thus, A^w = u)'(z) XAf{ z )i with xaa(^) defined as in (^6|), is the most general form of conformal 
variations defining normalizable zero modes. Since an arbitrary polynomial of degree Af depends on Af + 1 arbitrary 
complex parameters, the number of normalizable zero modes associated to conformal symmetry in the totally non 
degenerate case is 2Af + 2, just like in the totally degenerate case. 
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We conclude this section by observing that all the zero modes in the totally non degenerate case could have also 
been obtained through infinitesimal variations of the parameters defining u)(z). Under a variation of the parameters 
hi, the corresponding variation of w is A^lu = 5u/5bi = lu(z)/(z — bi). However, such a variation renders divergent I\ 
due to a non integrable singularity in z — hi, and infinitesimal variations of bi do not define normalizable zero modes. 
Under a variation of the remaining parameters luq and on we have 



Au uj = 5uj{z)/5uj = uj(z)/uj , 
A ai u = Su>/Sai — —lo(z)/(z — a-i) . 



(37) 



Substituting ( |37|) in (|28|), we have that both I\ and I2 are convergent, and eqs. (|37j) define 27V+ 2 normalizable zero 
modes. Obviously, these modes are not independent as they are precisely of the form (|35|), (fj6|). Actually, any linear 
combination of the zero modes associated to ([37|) can be used as a basis for the arbitrary functions xat(z) m (1^1) • 

Finally, also in this case one can easily check that the zero modes associated to ISO(2) symmetry as well as the 
gauge zero mode ( p6|) are contained in eqs. (^), (|36| ) (or, equivalently, in eqs. ([37])). 

Starting from these two limiting examples, it is easy to extract zero modes from any other intermediate choice of 
the arbitrary function uj(z) . In all the cases, the number of the normalizable zero modes is 27V" + 2 and it is thus only 
a function of the topological sector associated to the classical background. 

The thermal correction to the magnetic field due to the zero modes is given by 



b = en V dua v = -dd 



For instance, in the radially symmetric case, eq. (p8f) reads 



8W 2 



(ft) 



Me 



9nZ 



jAlo + wAu 
l + M 2 



— AT- 1 



(38) 



2Af 



1 + 



(ft) 



2Af 



(39) 



with 1— J\f < n < 1. It can be shown that the integral over the whole plane of eq. (|39j) vanishes for any 1—Af< n < 1 , 
so that zero modes corrections to the magnetic field do not modify the vorticity of the classical solutions, as expected 
on general ground. 

Having the zero modes corrections to the magnetic field, eq. (^) or (|39|), the solutions of the zero energy sector of 
the Schroedinger equation (E4T) straightforwardly follow. 



IV. CONFORMAL ANOMALY 

Although the overall factor 1/e 2 in eq. (|TJ) is clearly irrelevant at the classical level, one expects that statistically 
such a term is no longer inessential, as it introduces an explicit scale in the path integral. In this case the coupling 
constant will acquire a non trivial dependence on the arbitrary mass scale fi that any regularization precedure entails 
fl5f , and an "anomaly" is expected: the expectation value of the trace of the energy momentum is likely to be non-zero. 
Here we shall evaluate this expectation value by using the ^-function regularization technique. 

The determinant of the thermal fluctuation operator M, defined as the product of its non vanishing eigenvalues E(| 
Det M = Y[ n A n [M], is divergent due to the unbounded nature of its eigenvalues A„. A popular way to circumvent 
this problem, is to define the regularised determinant through the zeta function. The zeta function appropriate to 
our problem is defined by 

C(s\M)=Y l 'K S [M}= Tr (M- s ) . 

n 

The £ function is regular for Re s > d/m, where m is the order of the differential operator M and d the dimension of 
the manifold. However, £ can be analytically continued to the point s = and the regularised determinant is defined 

as 



" In Det ' M 1 > = lim 

J( = s^O 



-A C(S | M) 



= -C (0|M) . (40) 



Let us introduce the kernel 
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C (t;x,y|M) = 



r( s ) 



dbt 1 



/i(t;x,y|M)-P(°) (x,y) 



(41) 



where is the projector into the zero mode space and h ( t; x, y | M) the heat kernel satisfying the differential 
equation (d/dt + M x )h (t;x, y \ M) = with boundary condition h ( 0; x, y | M) = IS^ d \x — y). The heat kernel at 
x = y can be expanded in powers of t through the so called Seeley - de Witt |17]] coefficients 



h (t;x,x|M) = 



(4tt*) 



d/r. 



[a (x\M) +ai (x\M)t 



(42) 



Finally, standard manipulations |l^] permits to express the value of ((s = | M) in terms of the Seeley - de Witt 
coefficient a^/m and the number TV of zero modes of the operator M as 



C(s = Q|M) 



(47r)d/r. 



Tr / d d *a d/m (x)-N . 



(43) 



In our case, d/m — 1 and the total number of zero modes is N = 27V + 2, as seen in the previous section. 

It can be shown [fl9[ that the integral of the expectation value of the trace of the energy momentum tensor is given 

I'.v 



/ 



d 2 x(T pp ) =- 



Slogfj. 



(44) 



W = — lim Tr ' 

2 s^o 



[- 


fMV- 


ds 





[C(0) + In M 2 C(0)] 



where /i is the usual mass parameter introduced to render the operator (M/{j, 2 ) dimensionless, as required. Conse- 
quently, the integrated trace anomaly is just given by C(0), up to a sign, 



(T w ) = -C(0) . 



(45) 



In turn, from eq. (|43j), £(0) can be expressed in terms of the total number of zero modes of the operator M and in 
term of the integrated ai Seeley - de Witt coefficient. Thus, we only need to evaluate such a coefficient. 
If a (matrix valued) differential operator M is written in the form 



M = D+D^ 



X , x = x- 



(46) 



with D,, 



ldfj, — iCp, D+ = — ldfj, + iC+, C M = C+ , and X being matrices whose entries are solely classical 



fields (not operators), then the coefficient ai is just —X. In our case the operator M was already introduced in eq. 
(pp|), and after algebraic manipulations and taking the gauge d^A^ = and the self-duality conditions into account, 
it follows that M can be indeed written in the form (Eq), with 




(47) 



C-2 



i e 



( 01 2 



-0i A 2 

V-02 -A 2 



(48) 



X 



( 








Ln 


£12 \ 












L22 




ill 


£21 




-e 2 0i 02 


\ 


il2 


L22 


-e 2 0i 02 


-ie 2 2 - |e 2 2 ) 



(49) 



where 
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-eSabD^b 



(50) 



D^4> a — d^<j) a + eA^Eabipb being the U(l) covariant derivative of the matter fields <f> a . 
Consequently, the trace of a\ is 

tr a x (z) = -tr X = 2 e 2 + <#j) = 4 e J3(z) 



so that its integral is a topological invariant. From eqs. (45), (|43) it follows that 



(51) 



manifesting a trace anomaly for any Af ^ 1. 



V. DISCUSSION 



Among the local gauge theories admitting classical vortex solutions, the one we have presented is particularly 
interesting. First of all its classical solutions are very simple. All the solutions can be expressed in terms of an 
arbitrary analytic function, which is the arbitrary analytic function specifying the regular solutions of the Liouville 
equation. To any classical field configuration, there is an associated topological invariant (vorticity), which is the 
degree of the analytic function. The topological invariant is the same of that characterizing the solutions of the 2D 
euclidean NLSM, all of whose solutions can be also expressed in terms of an arbitrary analytic function. 

The trace of the energy momentum tensor vanishes on the classical solutions and, consequently, conformal trans- 
formations are a symmetry of the classical solutions: a conformal reparametrization relates two different solutions of 
the model. Related to this fact there is also the factorization property jli| ) of the coupling constant e, that is the only 
dimensional parameter of the model. Such a factorization is important also because it permits, for topological models, 
a temperature-like interpretation of the coupling constant ((3 — 1/e 2 ). In this context, the euclidean action should 
be interpreted as the potential energy of the system whereas the free energy should represent the energy fluctuations 
due to the interaction with the thermal bath (3. In the limit of small coupling (large (3), the partition function of the 
system is just the determinant of operator obtained byexpanding the classical action around the classical solutions up 
to the second order in the thermal fluctuations. The evaluation of the determinant requires a regularization procedure 
and, therefore, an arbitrary scale. Consequently, one expects that thermal fluctuations destroys scale (and conformal) 
symmetry, and in fact there is a conformal "anomaly" : the expectation value of the trace of the energy momentum 
tensor does not vanish, except for the special value Af = 1 (see eq. ([H])). When the classical background carries the 
minimum vorticity (Af = 1), conformal symmetry seems to survive, at least in the harmonic approximation. 

This is a very surprising property that certainly deserves a deeper analysis. At present, we do not know the exact 
reason of such a phenomenon. Nonetheless, there are indeed some properties that make the case Af = 1 different 
from the others. For example, the trace of the a\ Seeley - De Witt coefficient or, alternatively, the shape of the 
potential V — —eB cl felt by the Schroedinger particle (|24|), dramatically changes in the cases Af = 1 and Af ^ 1. Let 
us consider for convenience the radially symmetric case. Then, if Af = 1, the potential V — —eB cl = —e 2 (f>* cl (j) cl is 
a monotonic radial function, and therefore it has a single, non degenerate minimum at r = 0. On the contrary, if 
Af 7^ 1, the potential V has a maximum point at r = 0, whereas its minima are degenerate and located around the 
circle r = ro[(Af — l)/(Af + l)] 1 / 2Af . Thus, there is the intriguing possibility that the occurence of conformal anomaly 
could be related to the vacuum degeneracy of the free energy. 

Another feature that makes the case Af = 1 different from all the others is the following: the most general solution 
carrying vorticity 1 is obtained by choosing the arbitrary analytic function lo(z) — loq(z — a)/(z — b). It is easily 
recognized that such a function is in a one to one correspondence with an arbitrary transformation of the group 
SL(2, C), z — > £ = [az + b)/ (cz + d), ad — be = 1. This is a very special subgroup of the conformal group: it defines 
the projective transformations, that are the only conformal transformations providing invertible mappings of the whole 
complex plane onto itself. As a consequence of this fact, SX(2,C) is a kind of "residual" symmetry one has, once 
the the number of vortices and the total vorticity has been fixed: let the pair (Q,Af) denote an arbitrary background 
field configuration of Q vortices with total vorticity Af; then, applying an arbitrary SL(2, C) transformation, the 
background (Q,Af) is mapped into another background but with the same pair (Q,Af). By contrast, if a conformal 
transformation z — > not belonging to SL(2,C) is applied to a field configuration (Q,Af), it necessarily increases 
vorticity, i.e. (Q,Af) — > (Q' > Q,Af' > Af). Clearly, if Af = 1 (minimum vorticity) then necessarily Q = 1, and only 
in this case it happens that the most general element oj(z) characterizing the background (Q = l,Af = 1) belongs 
to the same group of transformations leaving the background (Q = l,Af = 1) unchanged. This observation is at the 
root to understand why the Af — 1 case is not anomalous. 
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Besides the conformal anomaly, we have also evaluated the zero modes of the thermal fluctuation operator. For a 
given fixed vorticity Af, we have derived the zero modes in two limiting cases of the classical background: the totally 
degenerate and the totally non degenerate cases, corresponding to a single vortex located at the origin and carrying 
vorticity Af, and to Af distinct vortices carrying each the minimum vorticity, respectively. All the other possible cases 
can be easily derived as intermediate between these two. Clearly, the explicit form of the zero modes dependson the 
classical background. However, the total number of normalizablc zero modes is only a function of the vorticity of the 
classical solutions. In addition, the contribution to the thermal fluctuations given by the zero modes does not change 
the vorticity of the classical solution. 

There are several aspects related to this model that deserve consideration for future investigations. An important 
problem is certainly a deeper understanding of the persistence of conformal symmetry in the Af = 1 case, and its 
possible relation with the non-degenerateness of the free energy vacuum. 

Concerning the eigenvalue problem (|2l]) or, equivalently, the determinant of the operator M, an exact evaluation 
seems very difficult. However, some approximate method should be available to investigate the system beyond the 
zero-mode sector, perhaps reducing the number of degrees of freedom by introducing some collective coordinates. 
Alternatively, one could try to investigate some particular limits that simplify the form of the matrix M . To this 
purpose, two particular cases should be mentioned: if ro — > 0, Liouville vortices becomes Aharonov-Bohm vortices, 
i.e. A cl ~ V0 and B cl ~ S(r). This limit greatly simplify the form of the fluctuation operator, although in this case 
the classical solutions become singular. Another interesting limit is the large Af limit; in this case the magnetic field 
has a significant non vanishing contribution only in a neighbourhood of r = r , and it could be easier to evaluate the 
determinant. 

Another interesting issue is the possibility that the model (|J) could be an 'effective action' of another, more 
elementary, model. This hypothesis could be supported by the fact that Liouville vortices always have an even 
vorticity in terms of the elementary flux quanta <&o = 27r/e when the Abelian coupling e is interpreted as electric 
constant. Consequently, it could be that the scalar field <fi is related to an electron pair condensate. A further point 
towards this direction is the one discussed at beginning of Section 3: in the very drastic approximation considered 
there, the partition function is the one of a Coulomb gas in its Debye phase that, in turn, can be equally described 
by a system of massive fermions. 

Finally, at the purely classical level, it could be interesting to investigate on the possibility of constructing non linear 
superpositions of Liouville vortices: just like the 't Hooft Polyakov monopole can be seen as a non linear superposition 
of infinite instantons equally separated in time pcfl , a non linear superposition of Liouville vortices could provide new 
soliton solutions of some lower dimensional gauge theory. 
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